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Abstract 

We explore the application and usefulness of the heavy quark symmetry to describe the weak decays of 
hadrons (mesons and baryons) containing two heavy quarks. Firstly we address the internal dynamics of a 
heavy-heavy bound system with the help of estimates based on potential models, showing an approximate 
spin symmetry in the preasymptotic quark mass region including charmonium, bottonium and B c meson 
states. However, no asymptotic spin symmetry should hold in the infinite quark mass limit in contrast 
to singly heavy hadrons. Predictions on semileptonic and two-body nonleptonic decays of B c mesons are 
shown. Furthermore, the stemming flavor and spin symmetries from the interaction between the heavy 
and light components in hadrons (combining in a "superfLavor" symmetry) permit their classification in 
meson-type supermultiplets containing singly heavy mesons together with doubly heavy baryons, and 
baryon-type supermultiplets containing singly heavy baryons together with some exotic doubly heavy 
multiquark states (diquonia). Exploiting their well-defined transformation properties under the superfla- 
vor symmetry group, we get predictions on the widths for some semileptonic and two-body nonleptonic 
decays of baryons containing both b and c quarks. 



To appear in Nuclear Physics B 



FTUV: 94-45 
IFIC: 94-40 



e-mail : mas@evalvx.ific.uv.es 
•Research partially supported by CICYT under contract AEN 93-0234 and by IVEI 



1 



1 Introduction 



History of Science teaches us that the discovery of new particles, far from complicating 
our world's description, opens new trends and offers a previously unknown insight into 
a more general understanding of the laws of Nature. In fact, a paradigmatic example in 
chemistry is the elements' periodic table whose full significance was not apparent until 
the development of the knowledge on the subatomic structure of matter and the birth of 
quantum mechanics. Analogously in elementary particle physics, the appearance of new 
(though shortly living) particle generations suggests or confirms new symmetries of matter 
and its fundamental interactions. Strange, charmed and beautiful hadrons have been good 
examples of it in the past and certainly still will be in the future. 

The existence of doubly heavy mesons and baryons (i.e. containing any pair of heavy 
quarks either c or b) can be viewed as a necessary consequence of the attractive strong 
interaction among quarks. So far, however, there has been no real possibility to produce 
them at an observable rate with existing accelerators, leaving aside charmonium and 
bottonium resonance families. 

Hopefully, this situation should definitely change with the advent of high-energy, high- 
luminosity hadron colliders like Tevatron and mainly LHC The two by now accepted 
experiments at LHC (ATLAS and CMS) include B-physics in their physical objectives 
and the feasibility of the detection of doubly heavy hadrons is being currently addressed 
|]| ||. Moreover, a third high-precision experiment at LHC dedicated to the study of 
Beauty particles could be finally approved. 

Recently the hadropro duct ion of heavy quarks in general, and of doubly heavy mesons 
and baryons in particular, has attracted a lot of theoretical work . Conversely, no such 
attention has been yet devoted to the analysis of the subsequent weak decay of these 
particles especially in the case of baryons, with the exceptions @| @ to our knowledge. 

At this stage, the purpose of this work is to make estimates on decay rates of doubly 
heavy hadrons in order to conduct an experimental search for their observation in a not so 
far future. We have focused, in particular, on initial state hadrons containing one b quark 
and one c quark. Indeed, in the case of mesons it represents the only possibility since the 
top quark would not form a bound state due to its large mass while for baryons the reason 
is an experimental one |J. We have tried to base our calculations on model- independent 
assumptions reducing whenever possible the input parameters within the framework of 
the approximations used in this work. Our final results, though presumably reasonable, 
ought to be viewed as a starting point for a more detailed treatment, as well as a checking 
limit of more exact derivations. 

This paper is organized as follows. In Section 2 we discuss the (partial) application 
of the heavy quark effective theory (HQET), primarily developed for heavy-light (h-1) 
systems, to the internal dynamics of bound states of two heavy quarks. In particular, 
the appearance of a spin symmetry at leading order (already pointed out elsewhere M- 



12] ) can be exploited in the analysis of the weak decays of the B c mesons. In Section 3 
we also study some weak decay channels of baryons with b and c quarks by means of a 
"superflavor" symmetry, obtaining predictions on decay widths and branching fractions 
whose numerical values are the first so far computed, to our knowledge. 
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2 Effective Lagrangians for heavy quarks 



In recent times, several attempts have appeared to describe the heavy quark strong dy- 
namics from more general grounds than particular models by simplifying the original 
QCD Lagrangian into a more manageable form. On the one hand, the non-relativistic 
QCD (NRQCD) is a systematic approach to the dynamics of heavy systems as increas- 
ing powers of the quark velocity inside the hadron [|L|. Once reformulated on a discrete 



space-time lattice, it constitutes nowadays an extremely active field of research fil . 

On the other hand, the last few years have witnessed a rapid progress in a new low- 
energy effective theory of QCD appropriate when dealing with hadrons consisting of one 
heavy component along with light degrees of freedom. Starting from the original La- 
grangian Cqcd = Q(itf — m o)Qi HQET can be formulated (T5| according to the following 
field redefinitions which remove the large part of the heavy quark momentum: 

h v ( x ) = e im ^ x P + Q{x) (1) 
H v (x) = J m Q v - x P_Q{x) (2) 

where P± = (l±^)/2 denote the usual projectors onto positive /negative energy states and 
h v and H v stand for the "large" and "small" components of the spinor field. In order to 
deal with hadrons containing a heavy antiquark, one must introduce the corresponding 



fields h„ and H„ 15 



The degrees of freedom corresponding to H v can be eliminated by means of the equa- 
tions of motion of QCD yielding a non-local effective Lagrangian. Then, a 1/ttiq power 
expansion in terms of local operators is achieved leading to 

= h^iv-Dhf + ^-[^ + Z(m Qi /fi)S^] + 0(l/ml) (3) 

where the superscript i labels the heavy flavor and those terms including 

tf» = h v {i >(iD ± )*h® (4) 
S (>) = f /^V^/^ (5) 

correspond to the (gauged) kinetic energy of the quark i in the hadron's rest frame, and the 
analog of the Pauli term describing the chromomagnetic interaction of the quark i with the 
light degrees of freedom. D± denotes the transverse covariant derivative, D± = D—v(v-D). 
The renormalization coefficient Z takes into account the short-distance physics in the 



chromomagnetic interaction [15]. 



The leading contribution in the Lagrangian of eq. (3) represents a static quark field, 
possessing two types of new symmetries not present in the original theory. In fact, the 
symmetry of this effective theory is larger than merely the product of the flavor and spin 
ones for the fermionic heavy quark fields, even extending to systems in which the massive 
colored particles have different spins. This issue will be discussed in more detail in Section 
2.2 and used in Section 3. 

2.1 Doubly heavy versus singly heavy hadronic systems 

In a weakly bound state of two heavy quarks such as charmonium or bottonium, the 
residual momentum k (defined as the quark's momentum measured in the hadron's rest 
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frame) neither should be too large on the average, amounting to a relative small off- 



shellness of the constituents [W\ 



Indeed, according to a large variety of potential models < k/rriQ > (directly related to 
the average internal quark velocity v' in non-relativistic motion) scales down with some 
power of the quark mass, i.e. < k/rriQ > ~ vtCq where v takes a value near 0.5 for a 



logarithmic-like potential, or near 0.34 for a "Cornell-type" potential fL7| . This behavior 
must be understood limited to a certain mass region becoming inappropiate in a dominant 
short-distance regime. 

Motivated by those arguments, we shall insist in applying the inverse heavy quark 
expansion generated by the field redefinitions (1) and (2) in order to get an effective non- 
relativistic Lagrangian to describe heavy-heavy (h-h) bound systems. As we shall see, 
important limitations of the low-energy effective theory will appear with regard to h-1 
systems, in particular the lack of flavor symmetry even at leading order whereas a spin 
symmetry will hold in a certain preasymptotic mass domain. 

In h-h systems the internal motion of the heavy constituents cannot be neglected even 
at leading order. This can be viewed as a consequence of the delicate balance between 
the average potential and kinetic energies in the bound state, as pointed out a long time 
ago in terms of the virial theorem [FS|. Thereby, one must keep the kinetic piece even at 



lowest order yielding a Schrodinger-type Lagrangian plus the Pauli term 

4{ = C ' ] + ^ Z ^l^ S0) + (VO (6) 

where 

C? ) = h^iv-D hf + ^- (7) 

The term including the operator (5) describes the chromomagnetic interaction of the quark 
% with the gluon field generated by itself and its heavy partner. 

Let us stress that in the Lagrangian expansion for a h-h system it is the order in the 
internal velocity v' rather than the dimension of the operators what matters in neglecting 
further terms Q Thus, from arguments based on power-counting rules (JOJ the expected 
magnitude of the Pauli term is suppressed by a power v' 2 with respect to C,. In fact, 
observe that for the sake of consistency in the characteristic magnitude, one should have 
considered some further terms of the expansion contributing at the same order in 
the internal quark velocity. Such next terms (up to order v' 4 ) would include relativistic 
corrections to the kinetic energy as well as spin-dependent effects like spin-orbit forces. The 
latter, however, vanish for s-wave states so that in this case the only relevant spin effect 
comes from the operator (5) while the former can be incorporated into a spin-independent 
Lagrangian beyond ■ 

Once taken into account, the chromomagnetic operator resolves the degeneracy between 
the spin-singlet and spin-triplet states, both in h-1 systems as in h-h systems. An important 
difference arises, however, in the behavior of the hyperfine splitting in the mg— >oo limit. 
In h-1 systems the degeneracy becomes asymptotically more and more exact. In sharp 
contrast, the validity of an approximate spin symmetry in doubly heavy states will be 



x This is the crucial difference with respect to h-1 systems where v' ~ Aqcd /rn,Q. 
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justified in a certain intermediate quark mass region, but it should certainly be broken 
above a more or less defined "threshold". We examine below this issue with the help of 
potential models. (Although potential models do not come directly from first principles 
they describe very well the general features of hadron spectroscopy). 

In order to assess the relative importance of the Pauli term in (6), we shall assume a 
Coulomb-plus-linear potential 

V(r) = V h or 

r 

For those states above the pure Coulomb well, the wavefunction at the origin is largely 
determined by the confining potential. In particular for the IS state, \ip(0) | 2 ~ e^m^a /2n) 
pOlj , where t\ = 1.935 and rhy = mQ/rriQ, / (tuq. + m^.) stands for the reduced mass of 
the QiQj system. For definiteness and simplicity, we ascribe the non-confining part of the 
potential to single-gluon exchange, hence k = cfol s where Cp = 4/3, 2/3 is a color factor 
for a quark- ant iquark pair, quark-quark pair respectively. We thus write the well-known 
expression [OT for the hyperfine mass splitting as 



Inserting the |^(0)| 2 expectation 



8cf ol s g 
3 m Qi + m Q . 



^M\ hyp * -r . . (9) 



and we see that AM\ hyp decreases whenever any rriQ i increases. In the particular case of 
J /if) or T, m Qt = m Qj = m Q and AM\ hyp ~ 16a s <j/9m Q . 

One can establish a relationship between the above expressions based on a potential 



model and an analogous definition of the A2 parameter for h-1 mesons [JloJ pl| , now for 
a doubly heavy bound system. To this end, we consider the spin-spin effect in the mass 
spectrum due to the chromomagnetic interaction of the quark i with a background of color 
field according to 

< H(v)\Jd 3 x 5C9) { (x)\H(v) > 

6M® = — ; r. m \ 

< H(v)\H(v) > 
where the correction term is here in particular 

and the hadron states are normalized non-relativistically. Let us now introduce 

<H(v)\S0\0)\H(v)> (0 

_ ( ) = d H A 2 (11J 

< H(v)\h * h$\H(v) > 

where dn = 3, —1 for the spin-singlet, spin-triplet state respectively. As far as we will 
be interested in s-wave states, the spin-spin interaction vanishes everywhere except in 
the origin x — 0. Thereby, taking into account that with the above normalization the 
denominator of (11) is unity, one can readily identify 6M®\ hyp = -d H Z\f /2m Qi . 
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In reality, the hyperfirie effect can be "ascribed" to either heavy quark i or j. Because 
of the symmetry of the spin interaction between both heavy quarks in a s-wave state 
(of the dipole-dipole type), the following identity can be formulated by requiring that 
SM(% yp = 6MW\ hyp . 

Z(m Q Jfii) mQ.Xi\ni) = Z(m Q J^j) m Qi Ar(/ij) = Z{m tj /^) myA 2 (/i) (12) 

where A 2 may be interpreted as resulting from the chromomagnetic interaction of a ficti- 
tious quark of mass rhij moving with the quarks' relative velocity through the background 
of the color field generated by both i and j quarks. Now, including a factor 1/2 to avoid 
double counting, we write in an explicitly symmetric form 

&M\ hyp = M H ,-M H = \[5M^\ hyp + 5M^\ hyp ] = 



Z{m Qi /iii) + Z(m Qj /nj) — A?Vi) = Z(rhij/n) A 2 (/x) 

m Qi rUQj J + m Qj 

(13) 

and thus ZX = (M^* — M|f)/4, in analogy to h-1 systems By comparison with eq. 
(9) we further get 

Zirhij/fx) A 2 (/i) ~ — ^(r^ 1 ) a (14) 

where r stands for the typical radius of the hadron, i.e. Tq 1 ~ cfol s {tq x )rhij. Note that 
A 2 is proportional to the string tension a and depends on the renormalization scale fi but 
Z(/j J )\2(fi) does not. Needless to say that eq. (14) must not be taken literally due its 
naive model dependence but rather as an indication of the close relationship between the 
A2 parameter and the string tension. Thus, the behavior of the hyperfine mass splitting for 
the J/ip and T, roughly of the type tuq 1 as in h-1 systems, amounts to the m^-independence 

of A 2 . 

A numerical estimate using a s ~ 0.4, a = 0.18 GeV 2 yields A 2 (/i^lGeV) ~ 0.13 GeV 2 , 
in significant accordance with the respective value of A 2 (/i— 2A) ~ 0.15 GeV 2 for a h-1 
meson. Indeed, this is not surprising for it could be directly guessed from the numerical 
similarity between the mass splittings of the D, D* states and the r] c , J/ip states . 

In summary, we advocate the validity of the application of the spin symmetry to the 
analysis of the weak decays of doubly heavy hadrons in analogy to singly heavy hadrons 
on the basis of the mg-independence and smallness of A 2 . 

Nevertheless, this emerging symmetry cannot be considered as asymptotically valid in 
h-h states. Since in the Coulomb regime v' ~ a s , we may then write AM\ hyp ~ niQa\. The 
"transition" mass domain where the linear behavior takes over the logarithmic fall off can 
be crudely estimated (from above) to be around m — e 4 AqcD — 10 GeV. Only below it 
(thereby including J/ip, B c and T states) spin symmetry can be used as an approximation. 
Indeed, in the language of potential models, a dominant Coulomb regime determines the 
wavefunction at the origin as |^(0) | 2 = (c^a^mjj) 3 / tt leading to the growth of the hyperfine 
splitting as follows from eq. (8). Stated in other way, A 2 would be proportional to the 
reduced mass squared yielding the growth of AM\ hyp . 
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2.2 Superflavor multiplets 

In the preceding Section we analyzed the interaction dynamics of a pair of heavy quarks 
borrowing to some extent the formalism developed for a h-1 system. In the following, we 
shall examine the interaction of a heavy diquark with a light component on the basis of 
an "additional" effective diquark theory. 

Let us start by recalling that QCD has, among others, a notable property: two color 
triplets can bind together to yield another color (anti)triplet. In particular, if both quarks 
are heavy enough, they will form an almost pointlike source of color behaving as a heavy 
antiquark with respect to a third light component in a hadron. This property of QCD, 
combined with the stemming spin-flavor symmetry at the very large quark mass limit 
should permit to classify hadrons containing one quark and two heavy antiquarks in mul- 
tiplets for the same velocity with similar properties as far as the light degrees of freedom 
are concerned. Indeed, the fundamental representation of the SU(6)<S> U(l) spin-flavor 
symmetry group consists of the two spin states of the fermionic quark field, the scalar 
di-antiquark field and the three spin states of the axial- vector di-antiquark field. All these 
objects appear to be the same to any bound light antiquark 0. In order to explicitly ex- 
hibit this symmetry, we shall write the representation of the superflavor group as "quark" 
supermultiplets (for each fixed v) expressed as nine-component column vectors 




* V 



where satisfies the constraint v^A^ = 0. (We have omitted any reference to the flavor 
of the constituents). The effective Lagrangian then reads 

= - H V M iv < D > tf„ 
where M. is a 9x9 mass matrix depending on the normalization of the fields which can 





be taken from J^] as 

M = 

and 

_ / 7° 

= oio 
V oos 

where g = diag(l, —1, —1, —1). 

One can form a superflavor "meson" multiplet by the tensor product of ^ v and one 
light antiquark field (there is no need to specify its flavor assuming unbroken light flavor 
SU(3)). Therefore, such hadronic supermultiplet puts together singly heavy mesons and 
doubly heavy antibaryons. The latter, though with different quantum numbers, really 
look like "mesons" in many respects, as commented more extensively in Section 3. In the 



2 



The underlying superflavor symmetry in the very heavy mass limit of the strong interaction was firstly developed by Georgi, 
Carone, Wise and Savage |22| 
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infinite quark mass limit, the behavior of the light degrees of freedom becomes independent 
of the heavy constituent flavors i, j and of the di-antiquark spin. Let us note that such 
supermultiplet is, however, not completely flavor-independent even at this asymptotic limit 
because of the mass dependence of the heavy di-antiquark internal dynamics 0. On the 
other hand, singly heavy baryons cannot be included in any meson supermultiplet since 
the groundstate light di-quark (either spin-0 or spin-1) should behave quite differently 
than a spin- 1/2 light antiquark. 

Instead, one can construct a baryon superflavor multiplet by combining either one 
heavy quark or one heavy di-antiquark with a light di-quark Moreover, we sort out two 
possibilities according to the spin of the light component. If it is an spinless object the su- 
permultiplet contains Aq baryons together with doubly heavy diquonia f\ (QiQjq q) . Other 
supermultiplets would include baryon (Eq, Eq) states with a spin-1 light component. 



3 Weak decays of hadrons with two heavy quarks 

Hadrons (mesons or baryons) containing two heavy quarks QiQj are particularly interest- 
ing since their lowest lying states would only decay weakly in contrast to charmonium and 
bottonium resonances. In the case of mesons, the only available possibility of gathering 
together two different heavy quarks is the B c meson. This is probably the reason why B c 
is nowadays object of growing attention both because of its weak decay and in hadron 
spectroscopy due to its interpolating position between the J/ip and the T p9f . 



3.1 Semileptonic decay of the B c meson 

As pointed out in the pioneering works by Nussinov and Wetzel |50j , Voloshin and Shifman 
|3lH the description of the decay of a h-1 hadron like a B meson simplifies dramatically 
in the infinite quark mass limit. When the heavy quark decays into another still massive 
quark, the surrounding brown-muck acts as a pure spectator with only one form factor 
(the Isgur-Wise universal function w = vi-v 2 ) governing the transition dynamics of 

the light degrees of freedom. At zero recoil C,(w = 1) is unity, reflecting a perfect overlap 



between the initial and final hadron wave functions [15 . 

On the other hand, it has been argued in recent literature about the possibility of 
testing the spectator behavior by experimentally discriminating among specific channels 
for b or c decays M in B c (be) mesons. Let us mention, however, that a simple-minded 
spectator picture of the B c decay should differ with regard to the analogous decays of 
B and D mesons. In h-1 mesons the spectator behavior makes sense as an asymptotic 



3 One may extend the often quoted pedagogical analogy between the atom and a h-1 meson. Now, the physics of the compound 
nucleus (the di-antiquark) is not at all the same for different isotopes although they have similar chemical properties (those of 
the light degrees of freedom). 

4 Henceforth, the reader himself will kindly realize whether we are referring to quarks or antiquarks so for the sake of brevity 
we generally suppress any explicit reference to antiquarks. 

5 Historically four-quark systems were in fact referred to as "baryonia" partly due to their expected decay into a baryon- 
antibaryon pair Q. In references p5| H Q the similarity between QiQjq q states and Aq baryons is established in the 
context of quark/potential models. Let us also mention that there is an alternative interpretation for the existence of such 
exotic states as di-mesons, i.e. (Qiq) — (Qjq) states p3]. Obviously, the latter would not be the supersymmetric partners of 
singly heavy baryons. 
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expectation whereas both heavy quarks are necessarily involved in the dynamics of a h-h 
hadron even for infinite constituent masses. 

Let us start by considering the b quark undergoing the decay, yielding the exclusive 
semileptonic channels: 

B c -> r)c{J/ip) +lu (15) 

Applying the spin symmetry discussed in Section 2 (see the appendix for more details), 
the expressions for the widths obtained in the non-recoil approximation are 

r(fl c -ifc + Zi/) ~ ^4 tp{x) f 2 l2 \v cb \ 2 (16) 



192 7T 3 



T(B C ^J/^ + tv) ~ -^_L [ V ( X ) + 0( X )] f( 2 \ Va \* (17) 

192 7T 

where x = m^jvai and </?, <fi are phase space factors for vanishing leptonic masses whose ex- 
pressions can be found in the appendix. Notice that we keep both </?, functions separately 
to facilitate later comparisons. The form factor is given by /x2 = \fm%~fm\ f]vikpi = V2), 
where the overlap factor f]\i is not fixed at zero recoil by a current conservation due to 
the lack of an underlying flavor symmetry. Instead, it can be estimated by means of a 
non-relativistic model of hadrons (see appendix) like the ISGW model according to 
expression (A.9). 

Next we want to apply our simplified approach to the decaying c-quark in the semilep- 
tonic channels of B c mesons: 

B c -> B S (B* S ) +lv (18) 

where the heavy quark spin symmetry and the non-recoil approximation will be employed 
once again. In || the authors made an estimate of the overlap integral between the 
BcS and _B s 's wavefunctions, denoted by Cl, by means of an operator product expansion 
retaining only the leading contribution. We have rederived the expression for Cl at zero 
recoil resulting, however, 

^4 (0 V2 / BsV ^:r o 3/2 (19) 

where r = l/cFrhb c a: s (mb c ) is the radius of the B c meson with rhb c denoting the reduced 
mass of the b and c quarks. Let us note some discrepancies between "our" formula and the 
original one presented in ||. On the one hand, there is a factor \/2 apparently due to a 
different normalization of states in (to 2v° in our case instead off ). Moreover, there is an 
extra factor \f?> in the denominator due to color f\ In fact, let us note that fs B y/mBs/^V^ 
can be interpreted as the wave function ipB s {fy at the origin according to quark models. 
Therefore, the "constant" ipB s {0) factorizes out in the overlap integral of the initial and 
final wavefunctions. (This can be accounted for since the variation of ipB s takes place 
over longer distances than the typical size tq of the B c ). However, let us observe that the 
evaluation of Cl in (19) should be somewhat overestimated because of the radial decrease 
of the B q wavefunction. 



6 In writing the initial B c state, a sum over colors must be understood in eq. (4.1) of ref. [g] and a normalization factor 1/V3 
should be included. Keeping color indices throughout, i.e. in the wavefunction and the color singlet weak current and taking 
into account the color Kronecker delta implicitly contained in eq. (4.4) of ||, the normalization factor l/\/3 survives all the 
steps of the calculation. 
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The decay widths come finally to be 



r(B e -+B. + lv) * ^4^)A 2 2i^l 2 (2°) 

19/ 7T 

r(B c ^B* s +l») ~ ^4^)+^)]/i2i^l 2 ( 21 ) 

192 7T° 

where / 12 = \Im?Jm\ Cl. 

Setting usual numerical values in (19) one gets from the above expressions widths of 
order 10~ 13 GeV, exceedingly larger than those expected from quark-model calculations, 
for example We therefore think that expression (19) actually overestimates the widths 
in the non-recoil approximation so we prefer to evaluate Cl employing the ISGW model 
according to eq. (A. 9) once again. 

In Table I we present our estimates turning out to be in remarkable agreement with the 
quark model predictions of ref. p3f for the same input parameters. (In fact, our results 
can be considered as oversimplified quark model calculations). This agreement extends 
as well over the results by other authors in spite of the simplicity of our approach. There 
is, however, an important discrepancy (common to all other references quoted in Table I) 
with the results by Bagan et al relative to the widths of those semileptonic channels with 
a final vector meson (either a J/f) or a B*). 



3.2 Two-body nonleptonic decay of the B c 

In this Section we address some two-body weak decays of the B c meson. Let us focus on 
the process induced by the decay of the 6-quark, 

B c -> Vc(J/ip) + vr (22) 

The exclusive channel corresponding to a final J / if) is one of the most promising from the 
point of view of its experimental detectability at a hadron collider 

We shall base our theoretical approach on two hypothesis. The assumption of the 
factorization permits to relate such decay with the semileptonic one at q 2 = m 2 ~ 0, 
according to 

r(fl c -7j c (J/^)+7r) 2f2 



dT/dq\ 2= JB c -^ Vc {J/^)+£v) 



6ttV; \VJC ~ 1 GeV< 



Besides, we keep the non-recoil approximation in the evaluation of the hadronic matrix 
element. Thus, making use of the set of equations (A. 6) and (A. 7) a straightforward 
derivation yields 

T(B c ^ Vc (J/iP) + tt) * G \ { ^ = ^ f 12 \V cb \ 2 (23) 

192 ir 6 mi 

where m 2 stands either for the mass of the r] c or the J /if). As seen from Table II, the 
agreement with other quark model predictions is remarkable. Furthermore, the general 
accordance manifested both in tables I and II suggests the extension of our simplified 
formalism to the weak decays of baryons with two heavy quarks. 
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3.3 Semileptonic decay of a doubly heavy baryon 

It has been pedagogically introduced in the context of HQET the analogy between an 
atom and a meson containing one heavy quark along with light degrees of freedom |34}] . 



On the other hand, in baryons with two heavy quarks (QiQjq), the massive quark pair 
should bind into a 3 source of color interacting with the light degrees of freedom as a 
static antiquark |37| . Assuming one-gluon exchange, the Bohr radius of the 

heavy diquark is f ~ 1/ 'c^a^r^ 1 )^,,- ~ 2r , where is the reduced mass of the (QiQj) 
subsystem. Indeed, one expects the size of the heavy diquark be presumably smaller than 
the typical length of the hadron, so that the analogous picture of a doubly heavy baryon 
in atomic physics would be a deuterium atom rather than an hydrogen molecule. 

If the two heavy quarks have different flavors there are two possible spin configurations 
of the compact diquark at lowest level. Supposing that the (QiQj) forms a spin-singlet 
(thus flavor- singlet) state the ground-state hadron will have spin 1/2 which we shall refer to 
as a Aij baryon. If the {QiQj) forms a spin-triplet (thus flavor-triplet) state the ground- 
state hadron can be either a spin-1/2 baryon (denoted by E^-) or a spin-3/2 baryon 
(denoted by E£ ). 

In order to apply the heavy quark spin symmetry to these doubly heavy baryons, we 
must realize two different scenarios: i) the "outer part" formed by the heavy diquark as 
a whole and the light component surrounding it, and ii) the heavy diquark subsystem 
itself. Ideally, unbroken heavy quark symmetries imply the mass degeneracy among all 
these hadronic states [^3| ||. Actually, the chromomagnetic hyperfine interaction should 
resolve the degeneracy leading to the mass level: 

miEij) < m(A^) < m(E*.) (24) 



according to . Let us remark that these authors apply the hyperfine splitting to doubly 



heavy baryons mimicking h-1 mesons, that is they only consider the spin interaction of the 
heavy diquark with the light component previously mentioned in i), obtaining 

3 1 

m(A^) = m(E i:j ) + -Am\ hyp = m(E*-) - -Am| %p 

Note that any discussion about the hyperfine splitting inside the heavy diquark has ap- 
parently been overlooked in and [23]. In fact, the asymptotic behavior of each chro- 



momagnetic interaction is completely different as stressed in the first part of this work: 
the former i) should fall off as any heavy quark mass increases whereas the latter ii) would 
finally violate the spin symmetry. 

Nevertheless, the inequality (24) should probably hold for c, b constituent quarks since 
the effect of the "inner" hyperfine splitting is suppressed by roughly a factor (1/2) 4 with 
respect to the "outer" hyperfine splitting (a (1/2) 3 comes from the the square of the 
(Coulombic) wavefunction at the origin and an additional 1/2 factor due to k = cf« s 
in equation (9)). Accordingly, observe that the lowest lying state seems likely to be the 
Eftc, i.e. the spin 1/2 baryon with a diquark which might be considered as a J p = 1 + 
quasi-particle. 

Finally, let us stress that among all the possible doubly heavy baryons to be produced 
at the new generation of hadron colliders, only those with one b and one c have a real 
chance to be observed. (This experimental issue will be discussed in more detail elsewhere 
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JT|). Hereafter, we shall concentrate on the weakly decaying spin-1/2 S fec baryon. 
Sftc — > N cc transitions 

According to our notation, if the heavy diquark (be) in the baryon undergoes a semileptonic 
weak decay yielding another heavy diquark (cc) (figure 1), the final groundstate baryon 
can be either S* c . Sometimes we shall design generically both of them as N cc . 

The quark-level current J M = C7 M (1 — 75)6 responsible for the hadronic transition 
T.bc^Ncc can be written as f|, 



J 11 = fj 12 (vvv 2 ) 



(25) 



where Aij(A\j) annihilates (creates) an axial- vector diquark (ij); 1712 (v\-v 2) stands for the 
wavefunction overlap between the parent and daughter heavy diquarks. 

Both initial Y, bc and final N cc particles belong to meson supermultiplets of the SU(6) <S> 
U(l) superflavor symmetry, as explained in Section 2. Thereby, their well-defined trans- 
formation properties under this symmetry group permits the use of the trace formalism 
ensuring a single form factor concerning the overlap of the light degrees of freedom at 
leading order approximation. We thus shall write the overall form factor 7712 governing the 
hadronic transition factorizing into two pieces || @, namely 

r)i2 = Vi2(vrv 2 ) €(vfV2) (26) 

At low recoil, the light degrees of freedom only play a marginal role in the decay, i.e. 
£(w ~ 1) ~ 1 as in the analogous transitions between h-1 mesons. Indeed, as shown in 
ref. 0, the light component does not change the overall decay rate of the diquark though 
channels the decay into different spin-states according to some selection rules, as a direct 
consequence of QCD in the limit of infinite quark masses. 

Since no flavor symmetry can be applied to h-h systems, neither r/ 12 nor consequently 
7712 are normalized to unity at zero recoil but can be calculated in a model-dependent way 
as in || using Coulomb wavefunctions, 

. , 1N 3/2 ( rh bc a s (rh bc ) m cc a s {m cc ) 

7712 {w = 1 = 8x 1 , [ x = r~ \ 1 1 - W2 

In the numerical evaluation we adopted m c = 1.5 GeV, m b = 4.8 GeV, a s {fh bc ) ~ 0.3 and 
Ois^Ticc) ~ 0.4. A more exact computation of 7)12 using wavefunctions from the Coulomb- 
plus-linear potential of reference |29| provides a very similar numerical value [], close to 
unity. 

Next, we apply the non- recoil approximation to the evaluation of rates for the semilep- 
tonic decays as explained in the appendix. One arrives at 

r[S bc ^S cc or Kc + * Mi) + 20(x)] A 2 2 \V cb \ 2 (27) 

57o n 6 

Discriminating S cc from E* c in the final state, one gets readily (see equations (17) and 
(18) in reference 0). 



r[E bc ^E cc + iu] 



G| ml 
576 7T 3 



—ip(x) + ~(f){x) 



/i 2 2 \V*\* {21 



7 I am indebted to P. Gonzalez for technical advice in this calculation. 
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r[E 



G% 



bcr 



2 

+ 3< 



fl2 \Vc 



eft 



(29) 



cc 1 576 7T 3 

In Table III we show the numerical results using the last expressions. It is worth to 
mention that the relative branching fraction of the semileptonic decays E;, c — >-E cc (E* c ) +£v 
is completely inverted with regard to the decay B c -^r\ c (J / 'ip) +lv . (The underlying diquark 
transition is different in each case, 1 — ► 1 in the former and — > 0(1) in the latter. The 
groundstate A cc is forbidden by the Pauli exclusion principle). 



J bc 



Nh transitions 



Let us assume now that the decay of the Eb c proceeds through the quark decay c^s, d 
yielding a final 6-flavored baryon. It might be conceivable that the diquark keeps up and 
the final baryon could be viewed as made up of a h-1 diquark and a light component |. 
Another interpretation assumes that the initial diquark actually breaks up: the b quark 
becomes a nearly static source of color whereas one light diquark is formed. We adopt 
this second (more standard) possibility henceforth supposing in addition that the Nj, lies 
at the groundstate. Still, the light dressing can be either a spinless or a spin-1 object. In 
the following we examine the former case, i.e. when the final baryon is of the A^-type. 

We are addressing here a transition connecting doubly heavy to singly heavy baryons. In 
order to take benefit from the compact trace formalism developed in |^2| for the evaluation 
of hadronic matrix elements, we need to recall the fact that both initial and final hadrons 
belong to (different) superflavor multiplets under the SU(6)(& U(l) spin-flavor symmetry 
group. Therefore, we shall use the following 9x4 matrix representations for the E;, c and 
At, (anti)baryons, 








\ 



V ul{v)Ca^v al5 ) 



u T A (v)C 





where C is the conjugation matrix and the spinors Ui are normalized to Ui(v)ui(v) 
satisfying fuiiv) = Ui(v). The corresponding bra matrices are defined as 



#( v ) = 7 °|rt( t 




2M 



(30) 



The trace formalism can be employed to compute the hadronic matrix elements for the 
Ef, c — >Nf, transition, as before. There are, however, some crucial differences with respect 
to the Efc— >N CC + £v decay as we shall see. 

The initial (be) diquark (with a very large mass ~ m& + m c ) should stand very close 
to the geometrical center of an spherically symmetric baryon. Furthermore, the 6-quark 
itself would tend to occupy a nearly central position inside the heavy diquark because of 
its relative larger mass. 

The mean distance of the b with respect to the center of mass of the (be) system (and 
to the center of the baryon as a whole too) is of order (rhb c /mb)f where f = 3/2a s rhb c 



This may make some sense in particular for an emerging s-quark which could be considered as semi-heavy to some extent; 
in fact, such a picture has been recently proposed as an alternative interpretation for any baryon containing a single heavy 
quark |3q]. 
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denotes the characteristic size of the (be) bound state, about twice the Bohr radius of the 
B c meson. Thereby, the typical distance of the 6-quark to the baryon center is ~ 3/2a s m b . 
On the other hand, the final heavy quark (i.e. the b) will stay near the central position of 
the final baryon as well, oscillating with typical "zitterbewegung" frequency 2m b and thus 
spreading over the corresponding length scale. Long-distance processes due to soft gluons 
give rise to residual momentum exchanges and intermediate quantum states with energy 
fluctuations of order k, leading also to a typical distance ~ l/kxk/m b = l/m b about 
the average position of the b quark. Observe that both of its initial and final oscillation 
amplitudes are of the same order (indeed they are l/m b effects) vanishing in the limit 
m b ^oo. 

Therefore, the decay of the c-quark at low-recoil occurs under the influence of a color 
field generated by a center-positioned, almost immovable 6-quark. (This contrasts with 
the decay Y, bc ^N cc where a reassembly of the final heavy diquark is necessary). Hence, no 
form factor relative to the heavy component will be considered in the hadronic transition 
amplitude, at least to leading order, or stated in other words, the wavefunction overlap 
corresponding to the 6-quark is assumed equal to unity. 

From the above discussion, one can heuristically deduce that solely those form factors 
related to the light degrees of freedom will come into play in the weak current. Accordingly 
the latter, J M = s~7 M (l — 75)0, acting on a S-type baryon can be written as 

J" = F^{v u v 2 ) h\A bcu 

F^ v is a tensor function depending on v\ and v 2 , which can be parametrized in general as 

F""(ui, v 2 ) = a(w) + b(w) « + c(w) v%v v 2 + d(w) ie^ v la v w 

where a, b, c, and d are universal functions concerning the light degrees of freedom. In 
particular at zero recoil (recall from the appendix the small range Vi-v 2 < 1.02 allowed in 
this process !) only the a-term gives contribution, all others vanishing since v"A bcu (vi) = 0. 
The a(u>~l) function absorbs the complexity of the non-perturbative decay dynamics 
of the light degrees of freedom near this kinematic point. Let us remark that it does 
not coincide with the Isgur-Wise function for transitions among members of meson-type 
supermultiplets. Nevertheless, a(w) is a true universal function no matters the flavor of 
the constituent heavy quarks. 

The hadronic matrix elements near zero recoil can be written as 

< A b (v)\r\i: bc (v) > ^ u A (v) 7^75 uj:(v) (31) 

Finally, we get in the non- recoil approximation v 1 = v 2 , 

T[^A b + Jv\ ~ [ip(x) + <P(x)] f 12 \V CS \ 2 (32) 



where f\ 2 = ym 2 /mi a(l). 

In order to have an estimate of a(l) one may conjecture that its numerical value should 
be intermediate between unity (i.e. the Isgur-Wise function at zero recoil) and the overlap 
of the light degrees of freedom corresponding to the D— > K transition (see figure 2). We 
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tentatively adopt the last wavef unction overlap factor obtained from the ISGW model 
with the parameters of Table A.I as an acceptable estimate for a(l). Indeed, note that 
the presence of the almost immovable fe-quark acting as an attractive center of color force 
should somehow shrink the size of the light degrees of freedom (especially compared to 
the Kaon system where such an static force center does not exist at all) thereby increasing 
their wavefunction overlap. However, this effect should be somewhat cancelled by the fact 
that the quark-quark interaction is about half as attractive as in the quark-antiquark case. 
Our predictions (with a large uncertainty) are shown in Table III. 



3.4 Absolute branching fractions 

In order to present our results in a more suitable way from the experimental viewpoint, 
we tentatively assume that the full width of the £f, c baryon is similar to the full width of 
the B c baryon. This assumption can be somehow justified by the fact that the A& baryon 



has similar lifetime as the B meson 40 



The main objection to this motivated analogy may arise from the fact the annihilation 
channel of the meson B c — > Iv has no analogous counterpart in the Af, c baryon decay. 
However, noting that the expected branching ratio of the B c into all leptonic channels is 
~ 18% |33| while the main part of the decay (82%) corresponds to "spectator" diagrams, 



it remains conceivable that the E^'s lifetime is not very different from the B c 's |33 
Therefore, adopting t(£{, c ) ~ 0.5 ps we get the set of branching fractions shown in Table III 
(since one would expect t(B c ) < r(Sb c ) we are underestimating these branching fractions 
in some sort). 



3.5 Two-body nonleptonic decay of the £ 



be 



In this Section we shall examine the two-body weak decay generated by the quark-level 
process b— >ccs 

S 6c -> E c + J/4) 

where E c denotes a ground-state Ag-type baryon. We have focused on this exclusive 
channel yielding a J/ip since its subsequent decay into a lepton pair could be useful to 
provide a tagging signal among the huge hadronic background at LHC |[|. Proceeding 
throughout in a very naive way, the hadronic matrix element 

< E c J/tp | J*% | S bc > (33) 

where J M = C7 M (1 — 75)6, = S7 M (1 — 75)0 can be Fierz-reordered. Thus assuming in 
addition the factorization hypothesis, it can be written as 

< ~ c I sY(l - 7s)& I S fec >< J 7V I c 7 M (l - 7s)c I > (34) 



We next relate this amplitude to the semileptonic decay at q = Trij,^ according to |}5[] as 
we did for the B r , 



dT/dq\ 2 _ m2 (S 6c ^ N c + tv) " ^ tj ^ SF lV - 



where sp stands for a suppression factor that we assume to be 4/9 (the factor 9 comes 
from color mismatch whereas the enhancing factor 4 is due to the existence of two final 
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c quarks to combine with the emerging c to yield the J/if>). Note that the denominator 
represents an "unreal" semileptonic process as an artefact with similar initial and final 
baryon states as the real physical process in the numerator. Then we use the following 
expression (see eq. (A. 7)) in the non-recoil approximation 



dT(l-^O) G 2 F [\ 3/2 + I2q 2 m 2 \ 1 ' 2 } 



fl2 \Vc 



,2 



cb\ 



dq 2 576 7r 3 m\ 

With regard to the overlap factors contained in / 12 we shall stretch the arguments used 
for the semileptonic decay £i> c — >N\, + tv, whereas this time the c-quark (always remaining 
within a region of dimension ~ l/^ c ) playing the spectator role instead of the b. There- 
fore, we (tentatively) disregard any form factor due to the rearrangement of the heavy 
component, whereas we evaluate the light component overlap according to the ISGW 
model like a ll B^K n transition. 

In table IV we present our predictions for the partial width and branching fraction of 
the process. It is interesting to note that the BF turns out to be of order 1%. Of course, 
such estimate suffers from many theoretical uncertainties but, however, is in remarkable 
accordance with an equivalent calculation by Mannel and Roberts [|I| for the analogous 
A&^A J /if) decay once taken into account the enhancing factor 4. 



4 Summary 

We have analyzed the dynamics of the weak decay of hadrons containing two heavy quarks 
employing the formalism of the Heavy Quark Effective Theory to some extent. Firstly 
we focused on the internal dynamics of heavy-heavy bound systems. Notice that con- 
trary to heavy-light systems, there is no flavor symmetry arising from the leading-order 
Schrodinger-type Lagrangian. 

We showed on the basis of estimates from potential models that in the range of masses 
limited to the preasymptotic domain including J /if), B c and T states, the Pauli chro- 
momagnetic term in the effective Lagrangian behaves as a genuine subleading term in a 
I/tuq expansion. For heavy quark masses larger than ~ 10 GeV, the Coulomb regime will 
prevail and the spin symmetry will be definitely lost. 

On the other hand, the interaction of a heavy diquark with a light component in 
a hadron exhibits a "superflavor" symmetry arising in the low-energy limit of the strong 
interactions, meaning that in the infinite mass limit both spin and flavor of the color source 
become irrelevant. Thereby, the massive quark and diquark fields can be joined into a 
common "quark" supermultiplet as far as their interactions with light degrees of freedom 
are concerned. Furthermore, such "quark" supermultiplet can be combined with either a 
light quark or a light diquark field to yield a classification scheme of hadrons according to 
the behavior of the light component. In particular, one can form a meson- type multiplet 
including singly heavy mesons and doubly heavy (anti)baryons, all of them with well- 
defined transformations properties under the superflavor group. All members of a hadron 
supermultiplet, though with different quantum numbers, would show similar spectroscopy 
concerning the excitations of the light component. In a parallel way, one can form two 
kinds of baryon-type multiplets by combination of the "quark" supermultiplet with either 
a spinless or a spin-1 light component. Thus, singly heavy baryons and exotic multiquark 
states (doubly heavy diquonia QiQjq q) would belong to the same supermultiplet. 
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In this work we made a phenomenological application of all the above ideas to the 
decays of hadrons containing b and c quarks in particular. Besides, we focused on those 
exclusive channels leading (at least) to a final charged lepton. In fact, to be realistic, the 
presence of a J/ip resonance decaying into a charged lepton pair will likely be an almost 
necessary condition to disentangle the decay from the huge hadronic background in hadron 
colliders (for more details see [[![]). 

Firstly, use was made of spin symmetry in the analysis of some weak decays of the B c 
meson. Predictions on partial widths and branching fractions by means of the non-recoil 
approximation are shown in Tables I, II. In general our evaluations are in agreement with 
other quark model calculations clearly favoring the final spin-triplet state versus the spin- 
singlet state. (However, a significant discrepancy to this general rule appears in the work 
of Bagan et al using QCD sum rules). 

Motivated by this general agreement, we have extended our formalism to the analysis 
of some weak decays of the £b c baryon. One ought to sort out two possible scenarios 
depending on whether the decaying quark is either the bottom or the charm quark. 

Let us suppose that the b quark undergoes the weak decay yielding a doubly heavy final 
baryon (iV cc ) and we look upon it as made of one heavy diquark along with light degrees 
of freedom. Then, both parent and daughter particles belong to meson-type multiplets 
and the trace formalism based on the superflavor symmetry developed by Georgi et al [P^] 
can be applied to straightforwardly. 

In this case, two form factors are required to describe the hadronic transition at leading 
order in the effective theory. One of them measures the overlap of the wavefunctions of 
the initial and final heavy diquarks, while the other is simply the Isgur-Wise function, 
common to transitions among members of this kind of supermultiplet. Moreover, near 
zero recoil the light component should not change the overall decay of the heavy diquark. 

Predictions on decay widths and branching fractions for semileptonic decays adopting 
the non-recoil approximation and assuming similar lifetimes for the baryon and the 
B c meson can be found in Table III. 

Instead, if the transition proceeds through the c-quark decay and the final baryon 
consists of the surviving 6-quark and one light diquark, we are actually confronted to a 
distinct physical situation. Indeed, the final hadron belongs to a baryon-type supermul- 
tiplet whereas the initial one is of the meson-type. We are thus dealing with a transition 
between members of different supermultiplets. 

Notice that although the heavy diquark would break up, the 6-quark would remain 
within a small central region of the baryon before and after the disintegration, in this 
particular channel. Therefore no form factor measuring the rearrangement of the heavy 
component in the baryon comes into play now at leading order. Instead, we introduced 
a form factor related to the light component in analogy to D^K transitions. Table III 
shows our predictions for the partial width and branching fraction. 

Lastly, we have evaluated the expected branching fraction for the two-body nonleptonic 
decay of the E^ baryon into S c and J/ijj relying on the factorization hypothesis. Our 
results are presented in Table IV. 
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Appendix 



In this paper, the hadronic matrix elements for weak transitions of doubly heavy hadrons 
are evaluated employing the non-recoil approximation (the parent and daughter hadrons 
are at rest in their common rest frame). We will retain, nevertheless, relativistic phase 
space factors in the calculation of decay widths but leptonic masses will be neglected. 

Throughout, subscripts 1 and 2 will design initial and final quantities. In particular 
mi, Vi, pi and rri2, V2, P2 stand for the mass, four- velocity and four-momentum of the 
initial and final hadron respectively. 

The validity of the non-recoil approximation as an estimate is strongly supported by 
the fact that the kinematic variable w = v\-v 2 is restricted to values close to unity. To 
appreciate this, let us consider a hadron containing two heavy quarks denoted by Qi and 
Qj. Let quark Qj, for example, undergo a weak decay into Qk- At q 2 = (pi — p 2 ) 2 = 0, 

(mi - m 2 ) 2 (m Qj - m Qk ) 2 

fr-**- = 1 + * 1 + 2 (mo , + mQ ,) (mQ , +mo J ^ 

Considering for definiteness the decay c— >s in an initial hadron containing one bottom 
quark/antiquark and one charm quark/antiquark, 

(bc) Vl -> (bs) V2 : (v r V 2 )rnax ^ 1-02 

for the mass values mj, = 4.8 GeV, m c = 1.5 GeV and m s = 0.3 GeV. 
If the b undergoes the weak decay, 

(bc) Vl -> {cc) V2 : (vyv 2 ) ma x - 1-29 

Analogous values should also hold for mesons and baryons with charm and bottom. Such 
small ranges of the velocities product should render the estimates of widths from hadronic 
amplitudes evaluated at w = 1 as an acceptable first approximation. 

On the other hand, the spin symmetry exhibited by eq. (7) permits to represent the 
covariant spin wavefunctions of doubly heavy bound states in a similar way as mesons 



with a single heavy quark, that is [15] 



J p = 0- : H{Q = -VM P+7v ; J V = oo~ : H(Q) = VM V4 (A.2) 
for a s-wave pseudoscalar meson. In the case of a s-wave bound state of two quarks, 

J p = 0+ : H(Q = VM V+ ; J v = oo + : H(Q = VM V+^ji (A.3) 

where P± = (l±^)/2 are energy projectors. 

In this way, the hadronic matrix elements corresponding to transitions among bound 
states containing two heavy quarks can be evaluated according to the well-known trace 



formalism: |39] 



< H 2 (v 2 )\Q k TQ J \H 1 (v 1 ) > = -r)i2{vyv 2 ) Tr^^TH^)} (A.4) 
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where the initial and final hadronic systems can be scalar, pseudoscalar, vector or pseu- 
dovector states; T is a combination of 7 matrices and 1712 plays the role of the dimensionless 
Isgur-Wise function for transitions between h-1 mesons |T5| , measuring as well the overlap 
of the initial and final hadronic wavef unctions. Yet, fji2(w) depends on the type of hadrons 
involved in the process and it is not a universal function of the velocity difference only. 
On the other hand, unbroken spin symmetry ensures a single 7)12 (reduced) form factor 
for hadronic transitions, irrespective of T, between different spin-states. Finally, let us 
stress that fju is not absolutely normalized to unity at zero recoil via the vector current 
since h-h systems do not belong to a multiplet of a heavy quark flavour symmetry group 
(there should be a mismatch between the parent and daughter hadron wavefunctions even 
at zero recoil). 

A summary of the relevant hadronic matrix elements in the non-recoil approximation 
(vi = t>2 = v ) for the weak process Qj^Qk W* is given by: 

< P , e 2 I U M I P > ~ ±2 7712 ^m~[m~ 2 

< l p , e 2 I A» I P > ~ ±2 77x2 

< P I A» I l p , e l > ~ ±2 r} 12 e? (A.5) 

< l p , e 2 I V I l p , e a > ~ =F 2 fju y/mm^ (e r e*) v» 

< l p , e 2 I A" I l p , e 1 > ~ =F 2 r)i 2 v/mT^ i e^ Q/3 ei Q e* 2/3 

where the upper(lower) sign applies when P = =p respectively; the vector and axial-vector 
currents are = Qkl^Qj an d A^ = Q^j^Qj. Upon contraction with the leptonic 
current ^7^(1 — 75)^, squaring, summing and averaging over polarization states whenever 
needed, the resulting differential widths for each massless lepton species and charge mode 
are ^ ^ 

-^(0 P Elector ~ ^ P 'Elector ~ ^3 /l2 l^f ( A ' 6 ) 

~ -— (o p -i p )i ~ ^(i p -i p )l ~ GUxV2 + 12 ^ A1/2 ) f \ V .A 

dtf K )lax * al -3dq2 [ )\ aX ial- 2dq 2^ ^ )\a X ial- 576ff 3 m 3 hi I V jk\ 

(A.7) 

where A=A(m 2 , m 2 , g 2 ) denotes the Kallen function, Vj^ is the KM mixing matrix element 
involved in the quark's decay. The (single) hadronic form factor / 12 at the non-recoil 
approximation is related to 7)12 by: 

fa = kI^- ma(w = 1) (A.8) 
V mi 

assumed nearly constant over the whole available q 2 range. This introduces an uncertainty 
lying within the framework of the approximations made in this work. Using the ISGW 



model p2[ at q max we may write 



^2 = M I 02 ( A - 9 ) 



Those values of the $ parameter of interest in our work are collected in table A.I. 
Integrating over q 2 between and q 2 max = (mi — m 2 ) 2 , 



max 

2 ™5 



no' -o p )L ctor ~ r(i p -.i p )|_ tor ~ ^ / 2 2 |u, fc | 2 (a.io) 
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^r(o^)|^ ^ r(i^^)|^ ^ ~r(i^)| asW - [¥>(*)+*(*)] A 2 2 |v iJk | 

(All) 

where ip and are phase space factors depending on the mass ratio x = m 2 /mi (neglecting 
lepton masses). Their expressions are: [|T3 

<p(x) = ~ [(l-x 4 )(l-8x 2 + x 4 )-24x 4 lnx] 

0(x) = 2 (1-a; 2 ) [(1 +x 2 ) 2 + 8a; 2 ] +24 (1 + a; 2 ) x 2 lnx 

In the limit x—*l, corresponding e.g. to the Shifman-Voloshin condition (mi — m 2 ) 2 << 
(mi +1712) 2 , i.e. small four-momentum transfer compared to the hadronic masses, phase 
space factors behave as: 

1 16 

-<P(x)~ v (x)^-(l-xf (A.12) 



In this extreme non-relativistic limit, 

r(o p -o^)|_ tor ~ r(i^i^)|_ tor * \ r{u p ^i p )\ axial ~ v{i p ^ p )\ axial * \ r(i p -i p ), a ^ 



G\ (mi -m 2 ) 5 2 . ,2 



(A.13) 



60 7T 3 

Assuming in addition that the spin-singlet and spin-triplet states are perfectly mass de- 
generate, the spin-rule yielding a factor 3 in favor of final vector states arises from the 
above expressions at once. In reality, however, broken spin symmetry may spoil such 
simple expectations mainly because of the different available phase space in each case. 
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Figure Captions: 



Figure 1) Graphs representing the weak decay of doubly heavy hadrons. For the 
weak interaction, the graphs may be interpreted as Feynman diagrams. 
In the case of the baryon the light quark should act as spectator at low re- 
coil, whereas the diquark can be considered as a decaying quasi-particle. 



Figure 2) A visualization of the hadronic transitions: a) £ 6c ^iV cc , b) £ 6c — >7V 6 , c) 
D^K. At zero recoil, the overlap factor of the light degrees of freedom 
in the case b) should be somehow intermediate between a) and c). 
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Table I. a). Partial widths of the semileptonic decays of B c (be) in units 10~ 15 GeV. The relative 
branching fraction refers to the final states with a vector meson with respect to those with a pseu- 
doscalar meson. The same input values of the parameters of the ISGW model, V c b = 0.046 etc as in 



ref. [33[] have been used to facilitate the comparison. 



Mode 


this work 


m 





H 


m 


B c ^> T] c + tv 


10 


10.6 


15 


14.2 


9.4 


B c ^ J/ip + tv 


42 


38.5 


44 


34.4 


12.2 


relative BF 


4.2 


3.6 


3 


2.4 


1.3 


B c ^ B s + tv 


18 


16.4 




26.6 


14.2 


B c ^ B* + tv 


43 


40.9 




44 


14.2 


relative BF 


2.4 


2.5 




1.7 


1 



Table I.b). Branching fractions from our results in Table La) using a lifetime of the B c equal to 0.5 
ps 
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Mode 


BF{%) 


B c ^ rj c + tv 


0.75 


B c ^ J/i/j + tv 


3.15 


B c ^ B s + tv 


1.35 


B c ^ B* + tv 


3.22 



Table II. a). Partial widths of nonleptonic decays of the B c in units 10 15 GeV, where we have used 
the same input values as in reference [Q. 



Mode 


this work 


PI 


m 


B c ^ r] c + n 


2.4 


2.1 


3.29 


B c ^ J /if) + 7T 


2.3 


2.2 


3.14 


relative BF 


~1 


~1 


~ 1 



Table II. b). Branching fractions from our results in Table II. a) using a lifetime of the B c equal to 
0.5ps[f53). 



Mode 


BF(%) 


B c —> Tlc + n 


0.18 


B c -> J /if) + 7T 


0.17 
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Table III. Partial widths of semileptonic decays of the T,j, c baryon (bcq) in units 10 
have used M Stc = 6.93 GeV, M Scc = 3.63 GeV, M s * c = 3.73 GeV, M Ab = 6.0 GeV |5 
estimates of the absolute BF we set the Ef, c lifetime equal to 0.5 ps. 



15 GeV. We 
For the 



Mode 


width 


BF(%) 




40 


3.0 


S 6c ^ £* c + tv 


12 


0.9 


S bc -^ A b + Iv 


34 


2.6 



Table IV). Width (in units of 10 15 GeV) and branching fraction of a two-body nonleptonic decay 
of the Sftc assuming its lifetime equal to 0.5 ps, M= c = 2.47 [[40| GeV and fj/^ = 0.385 GeV |33| . 



Mode 


width 


BF{%) 


S fec ^ E c + J/4> 


10 


0.8 



Table A.I. Values in GeV of the $ parameter in the ISGW model |32|] needed for the evaluation of 
several wavefunction overlaps at zero recoil used in this work. 



% = 


K 


D 


B 


B s 




B c 


Pi 


0.31 


0.39 


0.41 


0.51 


0.66 


0.82 
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